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ABSTRACT 

We study the importance of baryonic physics on predictions of the matter power spectrum as it is relevant 
for forthcoming weak lensing surveys. We quantify the impact of baryonic physics using a set of three cos- 
mological numerical simulations. Each simulation has the same initial density field, but models a different set 
of physical processes. The first simulation evolves the density field using gravity alone, the second includes 
non-radiative gasdynamics, and the third includes radiative heating and cooling of baryons, star formation, and 
supernova feedback. We find that baryonic processes alter predictions for the matter power spectrum signifi- 
cantly relative to models that include only gravitational interactions. Our results imply that future weak lensing 
experiments such as LSST and SNAP will be very sensitive to the poorly-understood physics governing the 
nonlinear evolution of the baryonic component of the universe. The net effect is significantly larger in the case 
of the model with cooling and star formation, in which case our results imply that contemporary surveys such as 
the CFHT Wide survey may also be sensitive to baryonic processes. In particular, this effect could be important 
for forecasts of the constraining power of future surveys if information from scales i > 1000 is included in the 
analysis. We find that deviations are caused primarily by the rearrangement of matter within individual dark 
matter halos relative to the gravity-only case, rather than a large-scale rearrangement of matter. Consequently, 
we propose a simple model, based on the phenomenological halo model of dark matter clustering, for baryonic 
effects that can be used to aid in the interpretation of forthcoming weak lensing data. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: clusters - clusters: formation - methods: 
numerical 



1. INTRODUCTION 

Contemporary determinations of the cosmic energy bud- 
get using a variety of cosmological probes all indicate that 
the majority of the energy in the universe (~ 70%) is in 
the form of dark energy with negative p ressure that drives 
an accelerated cosmic expansion (e.g., [Riess et al. 119981 : 
Perlmutter e"tai1ll999tjTeg mark et al. 2004; Riess et alJl2004t 
Eisenstein et alj l2005t ISpergel et alj 120061: iTegmark et alj 
20061: lAstier et alj|2006h IWood-Vasev et al.N2007l) . The un- 



known nature and properties of the dark energy is widely rec- 
ognized as one of the most important and fundamental prob- 
lems in cosmology, if not all of physics. 

Future large surveys, such as the Dark Energy Sur- 
vey 4 (PES), the S uperNova/Acceleration Probe 5 (SNAP, 
lAldering et alj 12004). and the Large Synoptic Survey Tele- 
scope 6 (LSST), are expected to measure the matter density 
fluctuation statistics from the linear t o the nonlinear regime 
with unprecedented precision (e.g., lAldering] 120051: iTvsonl 
120051) . The statistical power of such precision measurements 
should allow for s tringent constraints of the properties of 
dark energy (e g.. iHu & TegmarklfT999l: lHulll999t iHutererl 
20021: lHeavensll2003l: iRefregierl 120031: iRefregier et alj|2004t 
Song & Knoxl 120041: ITakada & Jain! 12004 ITakada & White! 
2004[ iDodelson & Zhangl 120051: lAlbrecht et all 120061: IZhanl 
2006). However, realizing this goal requires theoretical pre- 



dictions for the power spectrum in the linear and nonlinear 

1 Department of Astronomy & Astrophysics & Kavli Institute for Cosmo- 
logical Physics, The University of Chicago, Chicago, IL 60637 USA 

2 Enrico Fermi Institute, The University of Chicago, Chicago, IL 60637 
USA 

1 National Science Foundation Astronomy and Astrophysics Postdoctoral 
Fellow 

4 http : //www . darkenergysurvey . org/ 

5 http : //snap . Ibl . gov 

6 http: //www. Isst .org 



regimes and as a func tion of cosmology with an accuracy of 
order a percent (e.g., iHuterer & Takada 20051 iHuterer et alj 
120061) . 

The theoretical aspect of this program is daunting. To 
realize fully these goals, the power spectrum must be cali- 
brated on scales where nonlinear effects are important. Con- 
sequently, the calibration program relies on num erical simu- 
lation s of cosmological structure formation (e.g., Ann is et ail 
2005). However, it has not been yet demonstrated that mod- 
ern cosmological simulations can achieve the required accu- 
racy in the nonlinear regime. On the contrary, recent system- 
atic studies show differences of order ~ 10% between dif- 
ferent simulatio n codes at the rele vant range of scales (k x 
1 - 10 Mpc^ lHeitmann et al.l2005l) . For example, the limited 
numerical resolution of simulations, and therefore the limited 
ability of simulations to resolve the inner structures of dark 
matter halos and the properties of halo substructure, can give 
rise to non-n egligible systemati c effects on predicted power 
spectra (e.g., lHagan et al.ll2005l) . 

So far, systematic studies of the nonlinear power spectrum, 
including the calibration of popular analytic, phenomenolog- 
ical m odels, such as the halo model (see ICoorav & Shethl 
l2002l for a review), have utilized dissipati onless Af-body sim - 
ulations of gravitational clustering (e.g., ISmith et ai]|2003l) . 
However, in the observed Univers e approximately ~ 1 5% of 
matter is in the form of baryons (Spergel et al. 2006). Al- 
though baryons are sub-dominant and trace the distribution 
of dark matter at the onset of structure growth, the final dis- 
tribution of baryons in halos differs from that of dark matter 
significantly because they differ in their subsequent dynami- 
cal evolution. 

This differentiation has been observed routinely in cos- 
mological simulations of the formation of galaxy clus- 
ters. In general, these studies show that baryons have 
a more extended distribution than that of dark matter 
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(e.g., iFrenk et alii 19991) and this result is borne out by obser- 
vations (e.g., iDavid et alJH995b IVikhlinin et al.ll2006l) Re- 
cent simulations (e.g., iRasia et al.l2004tlLin et alT2 006) have 
shown that the transfer of energy from dark matter to baryons 
can modify the concentration of the dark matter radial dis- 
tribution by as 10%, even in simulations that do not include 
baryonic dissipation. Baryons may also be redistributed from 
high- to low-density regions via energetic AGN feedback, giv- 
ing rise to another, pote ntially non-negligible, effect on the 
matter power spectrum dLevine & Gnedinll2006l) . 

Furthermore, baryons dissipate energy through radiative 
processes, leading to the condensation of baryons in the cen- 
tral regions of halos where gas densities sufficient for star for- 
mation can be achieved. The condensing baryons pull dark 
matter along with t hem, leading to an increase in the dark 
matter density (e.g., IZeldovich et al.lll9 80; Blumenth al et al.l 
[1981 lGnedinetatll2004t ISellwood & McGaughll2005T) ~As 
a result, in any treatment of structure growth that self- 
consistently includes the condensation of baryons and the for- 
mation of galaxies, not only will some fraction of the total 
mass be redistributed into galaxies with comparably small 
spatial extent, but the dark matter halos themselves will have 
more compact, concentrated internal mass distributions. 

Previous analytic estimates of the effects o f baryon dis- 
sipation on the ma tter power spectrum (e.g., IWhitd l2~004t 
Zhan & Knox 2004) suggest that the effect should be confined 
to small scales. However, these preliminary studies are based 
on approximate, phenomenological models and may not have 
taken all of the relev ant effects into account. For instance, 
IZhan & Knoxl (|2004) model the relative influence of the hot 
baryons of the intracluster medium only, and neglect the in- 
fluence of a cold com ponent. 

A recent study by [Jing et al.l (120061) uses numerical sim- 
ulations to model self-consistently the interaction between 
baryons and dark matter and shows that dissipational physics 
can affect the power spectrum on scales k ~ 1 - 10/; Mpc" 1 
by > 5 - 10%. Our results provide qualitative confirmation 
of this result. Moreover, as we show below, in addition to 
the well-known increase of dark matter density in halo cen- 
ters in response t o baryon cond ensation, an effect included in 
the calculation of IWhitd (120041) . simulations indicate associ- 
ated changes in the distribution of matter at larger radii. This 
effect leads to enhanced concentration of the overall matter 
distribution in halos and affects the power spectrum on scales 
larger than previously thought. 

In the present paper, we use high-resolution, cosmologi- 
cal simulations of structure formation to study the effect of 
baryons and dissipation on the power spectrum of matter fluc- 
tuations at scales k ~ 0.1 — 10/iMpc" 1 . In particular, we per- 
form a series of three simulations each of which begins with 
the same set of initial conditions but includes different matter 
components and physical processes. The first is a collisionless 
Af-body simulation that models only a dark matter component 
and therefore includes only gravitational interactions. In the 
second simulation, we include a baryonic component evolv- 
ing according to an Eulerian hydrodynamics method. In this 
second simulation, the baryonic component is not permitted to 
cool radiatively. In the third simulation, we include radiative 
cooling for the baryonic component and prescriptions for star 
formation and feedback processes. We use the results of these 
simulations to study the relative impact of baryonic physics on 
the matter power spectrum and use the analytic halo model to 
show how these effects of baryons on the power spectrum can 



TABLE l 
List of Simulations 



Simulation 


Dark Matter 


Baryons 


Cooling/SF 


Resolution 


DM0 


yes 


no 


no 


0.9ft" 1 kpc 


DMG.NR 


yes 


yes 


no 


1.8ft" 1 kpc 


DMG.SF 


yes 


yes 


yes 


3.6ft" 1 kpc 



be understood in terms of differences in the density distribu- 
tions of baryons and dark matter in simulations with baryonic 
physics. We conclude with a brief discussion of the impli- 
cations of our results for upcoming precision weak lensing 
measurements of the matter power spectrum. For the purpose 
of inferring the properties of dark energy from weak lensing 
measurements, we suggest that it may be possible to encap- 
sulate the effects of galaxy formation in a small number of 
parameters. 

The paper is organized as follows. In §|2]we describe the set 
of cosmological simulations used in this study and in § [3] we 
describe the results of our analysis. These results include both 
three-dimensional and convergence power spectra as well as 
the structural properties of the dark matter halos and corre- 
sponding baryonic components in the simulations. We de- 
scribe the halo model for matter clustering in § [4] and use it 
to aid in the interpretation of the simulation results of § [3] 
In § |H we also suggest a method, based on the halo model, 
that can be used to model the baryonic physics present in our 
simulations. We summarize our conclusions and discuss their 
implications in §|5] 

2. SIMULATIONS AND ANALYSIS METHODS 

We investigate the influence of baryons on the matter power 
spectrum using cosmological numerical simulations. We sim- 
ulate the formation of structure in a cubic volume 60/t" 1 Mpc 
on a side in the "concordance" ACDM cosmology (Q.m = 0.3, 
Q A = 0.7, Q B /i 2 = 0.021, h = 0.7, tr 8 = 0.9). 

We use a set of three simulations, each of which begins with 
the same initial density field but includes different physical 
processes. The first simulation, "DM0," is purely dissipation- 
less and includes a collisionless dark matter component only. 
The second simulation, "DMGJJR," follows both dark matter 
and baryonic components but does not include radiative cool- 
ing for the latter. The gas in this run is thus modeled in the 
non-radiative (or "adiabatic") regime. The third simulation, 
"DMG_SF," treats the baryonic component including radiative 
cooling and heating, star formation, and feedback from super- 
novae. The inclusion of these physical processes allows for 
the formation of galaxies in the DMG_SF simulation as cooled 
gas forms a condensed component, a fraction of which is con- 
verted into stars. A basic description of these simulations is 
provided in Table Q] 

As the simulation volume is fairly small, cosmic variance 
and finite volume effects are significant at scales correspond- 
ing to k ~ 0.1 - lO/i Mpc -1 . This fundamental limitation 
prevents us from presenting absolute estimates of the power 
spectrum at precisions of order a percent. Additionally, we 
detect the effect of limited volume on the largest scales, 
k ~ 0.25/;" 1 Mpc, which is close to the scale at which density 
fluctuations have become nonlinear by z — (see Figure [T]). 
The transition to nonlinearity occurs between the fundamen- 
tal mode and this scale, so the finite box size may also af- 
fect the growth of the fundamental mode. However, these ef- 
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Fig. 1. — Upper Panel: The simulated power spectra for the TV-body 
(DM0), non-radiative (DMG.NR) and cooling (DMG.SF) simulations at z = 0, 
as well as the linear pow er spectrum obtained from the transfer function 
of lEisen stein & Hu 11999) and the actual initial density fluctuation spec- 
trum used in the simula tions. Note that the deficit of power relative to the 
lEisenstein & Hul (1999) spectrum is partially due to our initial realization, 
and partially due to nonlinear evolution. Lower Panel: The simulated power 
spectra now plotted relative to the initial density fluctuation spectrum. Note 
that this panel shows the power spectrum normalized to the initial power spec- 
trum realization used in our simulations and is not a residual compared to the 
Eisenstein & Hu 1 1999) spectrum. 

fects should be similar in all three simulations, and we are pri- 
marily interested in the relative effects of baryonic processes, 
which should be discernible through a comparison of the DM0, 
DMGJJR, and DMG_SF simulations. 

We perform the dissipationless simulation (DM0) us- 
ing the Adaptive Refinement Tree (ART) A^-body code 
( IKravtsov et all 1 19971: lKravtsov|[T999l) . The ART code em- 
ploys adaptive refinement in space and time in order to 
achieve the large dynamic range necessary to model the 
detailed structure of dark matter halos. In the DMGJJR 
and DMG_SF simulations, we follow the evolution of the 
gaseous baryonic component using an Eulerian hydrody- 
namics solver on the same adaptive mesh of the A^-body 
ART code, as described by Kra vtsov et alj d2002l) . How- 
ever, we perform the two simulations that include baryons 
utilizing the new, distributed-memory version of the N- 
body+gasdynamics ART code (Rudd & Kravtsov, in prepa- 
ration). 

The simulations with baryons treat both the collisionless 
gravitational dynamics of the dark matter and stars as well 
as hydrodynamical evolution of the gas component. Hydro- 
dynamic fluxes are computed at the boundary of mesh cells 
with a second order Godunov-type Riemann solver, which 
is particularly efficient at resolving shocks (typically within 
~l-2 grid cells). The simulation with radiative cooling treats 
several additional physical processes, including star forma- 
tion, metal enrichment and thermal feedback due to super- 
novae Type II and Type la, self-consistent advection of met- 
als, metallicity-dependent radiative cooling , and UV heating 
due to a cosmological ionizing background (Haardt & Madau 
1996). The cooling and heating rates take into account Comp- 
ton heating and cooling of the plasma, UV heating, and 
atomic and molecular cooling. The cooling rates are tabulated 
for the temperature range 10 < T < 10 9 K on a grid of metal - 
licities and UV intensities using the CLOUDY code (ver. 96b4; 



iFerland et aT]|1998l) . The CLOUDY cooling and heating rates 
take into account the metallicity of the gas, which is calcu- 
lated self-consistently during the simulation, so that the local 
cooling rates depend on the local gas metallicities. 

Star formation in the DMG_SF simulation is imple- 
mente d according to the observationally-motivated recipe 
(e.g., lKennicut u [l998h : p* oc p^/f,, with £» = 4 x 10 9 yrs. 
The gas in mesh cells is converted into collisionless stellar 
particles stochastically on timescales tsf = 10 s yrs, where 
the mass of the particle created from the gas is given by 
m» = p*TsFV ce ii, wnere Vceii is the volume of the cell. This 
prescription effectively averages the star formation rate in 
time, and is used to prevent large numbers of low-mass stel- 
lar particles from forming in high-resolution cells that contain 
little mass owing to their small volumes. The code also ac- 
counts for stellar feedback on the surrounding gas, including 
the injection of energy and heavy elements (metals) via stellar 
winds, supernovae, and secular mass loss. 

All three simulations follow the evolution of 256 3 dark 
matter particles (implying a particle mass of m p = 1.07 x 
10V M in DM0, and m p = 9.17 x lO 8 /^ 1 M in both DMGJJR 
and DMG_SF). Simulation DM0 uses a 256 3 uniform root grid 
and up to of eight refinement levels, giving a maximum co- 
moving spatial resolution of w 0.9/i~' kpc. The simulations 
DMGJJR and DMG.SF each use a 128 3 uniform root grid and 
allow for eight and seven levels of refinement, resulting in 
minimum comoving cell sizes of of 1 .8/i 1 kpc and 3.6/z~' kpc 
on a side, respectively. The simulations with baryons use 
somewhat coarser resolution compared to the DM0 run to com- 
pensate for the additional computational cost of these simula- 
tions, however the resolution of all three simulations is small 
compared to the scales we discuss. 

We identify halos in the simulations using a variant of 
the B ound Density Ma xima algorithm (BDM , iKIvpin et alj 
I1999T) . as described in IKravtsov etal] (I2004I) . We use the 
publicly-available smooth code 7 to compute the local den- 
sity at the position at each particle smoothed with a 24- 
particle smoothed particle hydrodynamics (SPH) kernel. We 
then identify halo centers with peaks in the smoothed dark 
matter density field. All particles within a search radius of 
rf = 50h~ l kpc of a peak are removed from further consider- 
ation as potential halo centers. The value of the parameter rf 
is set according to the sizes of the smallest objects we aim to 
identify in the simulations. After identifying halos according 
to the BDM algorithm, we calculate masses and other halo 
properties using all the matter components present in the sim- 
ulation (i.e., dark matter, gas, stars) within a virial radius, R wn . 
The virial radius is defined as the radius of the sphere, cen- 
tered on the highest-density particle in the halo, within which 
the mean density is a contrast A v i r (z) with respect to the mean 
matter density of the Universe. The value of A v h(z) is set 
according to the spherical tophat collapse model and can be 
computed quickly and a ccurately using the fitting formula of 
iBrvan & Normanl (fT998). In the concordance ACDM cosmol- 
ogy used in this work, A v i r (z = 0) ~ 337 and A v j r — > 178 for 
z > 1 . In what follows, we consider only distinct halos, or 
halos that do not lie within the virial radius of another, more 
massive halo. 

We measure the power s pectra in the simul ation vol- 
umes using the method of iJenkins et al.l (I1998L see also 
IKravtsov & Klvpinlll999i) . The power spectrum is computed 

7 http: //www-hpcc .astro .Washington. edu/tools/smooth. html 
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Fig. 2. — The simulated power spectra for the non-radiative (DMGJJR) 
and cooling+star formation (DMG.SF) simulations are plotted in left and right 
columns respectively. In each case we plot the fractional difference of the 
power spectra from that of the dissipationless (DMO) simulation. We present 
results for a range of redshifts z = 0, 0.5, and 1.0, as labeled in the left panel 
of each row. In each panel we show results for all matter (solid line) as well 
as for the dark matter (dotted line) and gas (dashed line) components alone. 
Stellar mass in the (DMG.SF) simulation is highly biased relative to the total 
matter distribution and is omitted for clarity. 

in a series of wavenumber ranges [k u k 1+ \\ as follows. We first 
divide the computational volume into i 3 cubic subvolumes. 
We superpose the particle distributions in each subvolume 
and assign particle densities to points on a 512 3 grid using 
this composite particle dist ribution according to the clo ud-in- 
cell interpolation scheme dHocknev & Eastwood 1981). For 
the gaseous component, we assign densities directly from the 
adaptive refinement mesh. The fast Fourier transform of the 
density field then yields an estimate of the power in modes 
that are periodic with period L\, ox /i, where Lbox = 60h~ l Mpc 
is the comoving box size. The wavenumbers at the boundaries 
of each interval are set to k\ - 512/(27r/L box )/18 in order to 
minimize the difference between adjacent segments. We use 
i = 2'" with m = ... 4. We make no corrections for parti- 
cle shot noise, which is negligible on all of the scales that we 
consider. 

3. SIMULATION RESULTS 
3.1. Power Spectra 
3.1.1. Three-Dimensional Matter Power Spectra 

Figure|2]shows the difference between the simulated power 
spectra in the radiative (DMG_SF) and non-radiative (DMGJJR) 
hydrodynamical simulations relative to the spectrum mea- 
sured in the A^-body simulation (DMO). The spectra are plotted 
from the fundamental mode at k — 2n/L\, 0X as Q.lh Mpc -1 to 
k w 30h Mpc -1 below which the power should not be affected 
directly by numerical resolution. 

First, consider the power spectra from simulation DMGJJR 
with non-radiative gas physics. Relative to the A^-body sim- 
ulation DMO, the total matter power spectrum in the DMGJJR 
simulation is enhanced by of order ~ 5 - 10% at wavenum- 
bers 1 < £/(/iMpc~') < 10 at all redshifts. Decomposing 
this change into the contributions from the gas and dark mat- 
ter components, we see that the dark matter power is enhanced 
by ~ 10-15% on these scales. Conversely, the gas component 
exhibits diminished power relative to the A^-body simulation 



on all scales k > 0.7 h Mpc 1 . In the range of wavenumbers 
1 < fc/(/zMpc ) < 6, the power in the gas component is re- 
duced by roughly ~ 10% at all redshifts, but the power drops 
precipitously for k > 6/zMpc . The net result is that the 
total matter power spectrum in the non-radiative gas simula- 
tion is enhanced on scales 1 < k/(h Mpc -1 ) < 10, but drops 
rapidly for larger wavenumbers. These results are in qual- 
itative agreement wi th the resu l ts of the non-radiative sim- 
ulation presented by IJing et all d2006l) . but t here are differ- 
ences in detail. In particular. I Jin g et all (l2006h report that the 
gas power spectrum is reduced by ~ 20% or more on scales 
1 < fe/(/iMpc -1 ) < 6 at z = 0. Additionally. Uing et all (120061) 
report a slightly larger enhancement in the dark matter power 
on these scales. However, the net result is a total matter power 
spectrum that differs from their A^-body results by less than 
~ 5% for k < 6/iMpc -1 . While these two results are similar, 
it is important to note that they are achieved via different com- 
binations of suppressed gas power and enhanced dark matter 
power. 

Relative to the non-radiative simulation, the effect on the 
power spectrum in the simulation with gas cooling and star 
formation (DMG_SF) is significantly more dramatic. In this 
case, the gas power spectrum is greatly reduced on large 
scales (k < 10/iMpc -1 ) relative to the A^-body matter power 
spectrum. This large-scale bias of the gas is due to the 
fact that smaller, and therefore relatively more weakly clus- 
tered, halos have larger gas fractions then their larger and 
more strongly-clustered counterparts (see § 13.21 1. At k > 
10/iMpc~', the gas power spectrum rises dramatically, indi- 
cating that gas has cooled and condensed into dense clumps 
of cold gas within dark matter halo centers. Both the dark 
matter and total matter power spectra rise dramatically on 
scales k > 1/iMpc -1 . These scales approximately corre- 
spond to the sizes of the largest halos in our simulations 
(M ~ a few xlO 14 hr l M ) and this dramatic feature re- 
flects the fact that gas has cooled and condensed in halo cen- 
ters causing a conco mitant contraction of the dark matter ha- 
los themselves (e.g.. IZeldovich et alj|19 80; Blum enthal et alj 
fl98H lGnedm^tal]l200llSellwood & McGaughll2005l) . We 



elaborate on this point in § [ 

Notice that the total matter power spectrum lies above both 
the gas and dark matter power spectra. The total matter power 
spectrum includes a stellar component not depicted in Fig- 
ure [2] which is significantly more clustered on these scales 
than either the gas or the dark matter. Again, the effects of 
dissipation on the matter powe r spectra that we find agree 
qualitatively with the results of IJing et alj (|2006). However, 
we find a larger relative enhancement in both the dark mat- 
ter and tota l matter power sp ectra in the cooling case than is 
reported bv lJing et alj d2006l) . Also, we find a sharp increase 
in the gas power spectrum at k > 10/iMpc -1 in the DMG_SF 
si mulation which is not seen in the corresponding simulation 
of Jing et al](l2006l) . The source of these differences is not en- 
tirely clear, but could be attributable to the specifics of the im- 
plementations of cooling, star formation, and feedback in the 
numerical codes used, to different gaseous equations of state 
or, at least in part, to sample variance. Finally, we note that the 
cooling simulation exhibits a bias of ~ 2 - 5 % in the total mat- 
ter power spectrum at the largest scales (k < Q.3h Mpc ). We 
are unable to explain this bias quantitatively and it is observed 
neither in our no n-radiative simulation nor the simulation of 
IJing et all (2006). Finite volume along with the different non- 
linear evolution of the simulation with baryonic condensation 
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are likely culprits for the small differences in clustering on 
large scales. 

3.1.2. Convergence Power Spectra 

While it is simplest to examine changes to the three- 
dimensional power spectrum in the simulations, weak lensing 
surveys will measure the convergence power spectrum, which 
is a projection of P(k) convolved with a lensing weight func- 
tion. Assuming a flat cosmological model a nd using Limber's 
approximation (Lim b"eriri954t iKaiserll 19921) . the convergence 
power spectrum is given by 



n({+ 1) 



H(z) /A ' 



(1) 



where A 2 (k, z) = k 3 P(k, z)/2n 2 , the lensing weight is W(z) = 
3Q.MH 2 l g(x)(l + z)/2, andx is the comoving distance to red- 
shift z. The function g(x) = X^ d X' n (x)(x' ~X)/x'< where 
n(x) describes the distribution of lensed source galaxies nor- 
malized so the integral over all redshifts is unity, J n(x) dx = 
1 . The distribution of source galaxies is typically defined per 
redshift interval, and we adopt the notation n(z) as no confu- 
sion should arise. 

We illustrate the importance of the effects of baryonic 
physics on the convergence power spectrum as follows. We 
assume for simplicity that the source distribution is given by 
a thin sheet of sources at z s = 1, so that n(z) - S(z - z s ). 
Denoting Xs( z s) as me comoving distance to redshift z = 
z s = 1, this assumption gives g(x) = x(Xs ~ x)/Xs- Us- 
ing more common source galaxy redshift distribu tions, such 



as n{z) oc z exp(- 



/z Q ] 2 ) as in M a et alj d2006l ) or n(z) oc 



z 2 exp(-[z/zo]) as in iHut erer (20021 zo ~ 1 in both cases), 
affects the relative convergence power spectra that we quote 
at a level that is small compared to the baryonic processes 
that are the focus of this paper. We adopt a thin sheet of 
sources for simplicity as this choice makes our relative con- 
vergence spectra directly comparable to those presented by 
Jing et al. (2006), yet it remains a good representation of the 
relative spectra with more common source galaxy redshift dis- 
tributions. We compute the integral in Eq. (HJ by tabulating 
A 2 (k, z) on a grid in the space of wavenumber and redshift us- 
ing the outputs of the three simulations from z = to z = 1. 
There are 41 outputs available from the A^-body simulation, 
30 outputs available from the non-radiative gas simulation, 
and 34 outputs from the simulation with cooling and star for- 
mation. We extrapolate outside of the range of wavenumbers 
probed by our simulations using the halo model, the details of 
which are described in §|4] 

We compare this result to the statistical uncertainty in the 
measurement of the convergence power spectrum for both 
contemporary and long-term, future weak lensing surveys. 
We include both sample variance and intrinsic shap e noise 
under the assumption of a Gaussian density field (e.g., Kaiser 
1998), 
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where / s k y is the fraction of sky coverage of the survey, h is 
the effective number density of source galaxies on the sky, 
and y 2 is the rms intrinsic galaxy ellipticity. Note that the 
error assumes Gaussian statistics whereas non-Gaussian ef- 
fects are likely to be considerable over a large part of the 
range of scales where baryonic effects are also important 
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Fig. 3. — Fractional difference of the convergence power spectra derived 
from the non-radiative DMGJJR (solid line) and cooling DMG.SF (dashed line) 
simulations with gasdynamics relative to the power spectrum of dissipation- 
less iV-body simulation DM0. The three shaded bands correspond to the sta- 
tistical eiTors expected from the CFHT Wide Survey, SNAP, and LSST (from 
outermost to innermost band, respectively). Note that these error bands are 
based on the assumption of Gaussian statistics and do not include systematic 
uncertainties. Scales left of the dotted vertical line at C = 3000 are those 
scales typically used to forecast constraints on dark energy parameters from 
future weak lensing surveys. This demarcation is meant only as an approxi- 
mate guideline. 



(ICoorav & HulHOQlt ISemboloni et al. 2006}. As a near-term 
observational program, we consider a survey like the Canada- 
France-Hawaii-Telescope (CFHT) Legacy Wide Survey 8 cov- 
ering a fraction / s k y =4x1 0~ 3 of the sky, with a number den- 
sity of source galaxies of h — 13 arcmin~ 2 . We also consider 
the relative importance of baryonic physics compared to fu- 
ture weak lensing surveys such as LSST and SNAP. For con- 
creteness, we take f^y = 0.5 and h — 50 arcmin~ 2 for our fidu- 
cial LSST-like survey and f^ y — 0.025 and h = 100 arcmin~ 2 
for our fiducial SNAP-like survey. We adopt y 2 = 0.22 for 
the variance in intrinsic source galaxy ellipticity and a band- 
width of Alji — 1/10. We emphasize that these are estimates 
of statistical uncertainties under the assumption of Gaussian 
statistics and include neither deviations from Gaussianity nor 
systematic uncertainties. 

The resulting fractional difference of the convergence 
power spectra from that of the dissipationless A^-body sim- 
ulation are shown in Figure [3] The relative difference in 
the convergence power spectrum in the cooling simulation is 
as dramatic as was seen for P(k). For both the contempo- 
rary CFHT Wide Survey and the future LSST- and SNAP- 
like surveys, the systematic difference is far greater than sta- 
tistical uncertainties for £ > 10 3 . In practice, dark energy 
constraints are derived using multipoles less than some € max 
where deviatio ns from Gaussiani t y are assumed to be unim- 
portant (e.g., IWhite & Hull2000t IVale & White! 120031) . The 
vertical line at I — 3000 in Fig. [3] represents a typical value 
of Anax and is meant to serve as a rough guideline for scales 
relevant to forecasts that appear in the literature; ho wever, a 
wide range of choices for l max have been made (e.g.. iHutereri 
12001 iRefregieri 120031: Hefregier et alJ 12001 iMa et al.l " l2006h 
Izhanll2006l) . 

One might argue that current hydrodynamical simula- 
tions are not yet up to the task of predicting accurate 

8 http : //www . cfht . hawaii . edu/Science/CFHLS/ 
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Fig. 4. — Upper Panel: Cumulative mass function of dark matter halos. We 
plot the number of halos with virial mass greater than M as a function of M. 
Lower Panel: The fractional deviation from the dissipationless mass function. 
We plot the difference between the number of halos in the hydrodynamical 
simulations and the N-body simulation normalized by the number of halos in 
the W-body simulation, [N(> M) - Nmo(> M)]/N m0 (> M). 

power spectra due to notable shortcomings in modeling 
galaxy formatio n, such as the well-do c umented "overcooling" 
problem (e.g.. [Katz & White! Il993 | : [Suginohara & Ostriker 
|1998tlLewis et aUlgOOOt iPearce et alJl2000t iDave et alJl2001b 
Balogh et al. 2001; Borgani et al. 2002). However, notice that 
even in the case of the non-radiative simulation the power 
spectrum deviates from the A^-body case at a level that will 
be important to future efforts such as SNAP and LSST. This 
result suggests that our inability to make robust predictions 
for the evolution of the baryonic component of the Universe 
introduces a large, systematic uncertainty into the theoretical 
predictions for the growth of perturbations, and ultimately for 
the analysis and interpretation of such data. 

3.2. Simulated Halo Properties 

Although density perturbations are small and grow accord- 
ing to linear theory on large scales, on the small scales probed 
by weak lensing surveys typical density fluctuations approach 
or exceed unity and fluctuations on these scales evolve non- 
linearly. In fact, the power spectrum at k > 0.5/iMpc 1 
is dominat ed by the structures of v i rialized dark mat t er ha- 
los (e.g.. IScherrer & Bertschingerl fl99lt iMa & Frvl 120001 
IZhan & Knoxll2004l) . Given that the differences in the power 
spectra discussed in § 13.11 are primarily confined to these 
scales, we find it instructive to explore changes within these 
highly-nonlinear, bound objects. In this section we present 
the systematic differences in halo properties between the DM0, 
DMGJJR, and DMG_SF simulations. 

3.2.1. Halo Mass Functions and Gas Fractions 

We begin with the abundance of halos as a function of halo 
mass. Figure |4] shows the cumulative mass functions for the 
three simulations. The mass functions in the A^-body simu- 
lation and the non-radiative simulation are very similar for 
masses above ~ 10 12 h M . The two mass functions differ 
at lower masses due to differences in numerical resolution of 
the DM0 and DMG_NR runs. Consequently, we only consider ha- 
los with virial masses greater than 10 12 h~ l M Q (« 1000 dark 
matter particles) when comparing halo properties. 



Fig. 5. — The mean fraction of mass in baryons within the virial radius 
as a function of halo mass for halos in both the non-radiative simulation 
{dash-dotted curve) and the cooling+star formation simulation (solid curve). 
The baryon fraction is expressed in units of the universal value assumed in 
our simulations Qb /S^m- Note that the halo baryon fractions in the non- 
radiative simulation are below the universal value (exhibiting the so-called 
baryon bias), while the baryon fractions in the simulation with cooling and 
star formation are above the universal value. Hashed regions correspond to 
the 68% intrinsic scatter amongst halos, while error bars correspond to the 
estimated error of the mean. 

The cumulative mass function in the DMG_SF simulation is 
consistently ~ 10% higher than in the DM0 simulation over a 
wide range in mass. The strong redistribution of mass in the 
DMG_SF halos (see § 13.2.2b causes the halo mass (defined at a 
fixed overdensity) to be larger for the same halo between the 
DM0 and DMG_SF simulations, effectively shifting the cumula- 
tive mass function to larger masses. The clustering of these 
halos remains unchanged, however, which leads to a signifi- 
cant change in the halo bias-mass relation, a rapidly varying 
function of mass on cluster scales. 

In the two hydrodynamic simulations (DMGJJR and DMG_SF), 
we can also track the manner in which baryons are appor- 
tioned to dark matter halos. Figure [5] shows the baryon frac- 
tion as a function of halo mass normalized to the universal 
baryon fraction (Qb/^m = 0.143 in our cosmology). The 
halos in the non-radiative simulation exhibit baryon fractions 
that are « 5 — 10% below the universal value, consistent 
with previous simulations without cooling (|Eke et alJ 1 19981 : 
Frenk et alJ[l999t lEttori et~all 120061: iGottloeber et alJboQfj 
Crain et afll2006l) . Infalling baryons are heated by shocks 
at distances of order the virial radius or larger, exchange en- 
ergy with the dark matter component during relaxation, and 
thus no longer follow the same distribution as the dark mat- 
ter. Instead, the gas follows a more diffuse, extended pro- 
file and never achieves extremely high densities in halo cen- 
ters. Rather, as we discuss below, the gas profiles in the 
DMGJJR simulation exhibit cores of nearly constant density 
that extend to large fractions of halo virial radii. Note that 
the slight decline in the halo mass-baryon fraction relation to- 
ward M v h- ~ 10 12 h~ l M & may be due to limited numerical 
resolution. We leave to future studies the question of whether 
the mass-baryon fraction relation truly varies with mass in the 
non-radiative regime or whether there exists a fixed baryon- 
bias which is constant with mass. 

Halos in the DMG_SF simulation, on the other hand, have a 
~ 5 — 10% excess of baryons relative to the universal value. 
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This excess is caused by the ability of gas to dissipate kinetic 
energy and condense at the center. These baryons are then 
unlikely to make large excursions away from the halo center 
and they make the potential well seen by newly accreted ma- 
terial considerably deeper than in the A^-body or non-radiative 
cases. 

3.2.2. Halo Density Profiles 

Halos in dissipationless A^-body simulations have 
spherically-a veraged density prof iles that are well described 
the profile of iNavarro et all (1 19971 NFW), 
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The relative concentration of mass toward the halo center is 
described by the concentration parameter c. The concentra- 
tion parameter has been studied by numerous authors and the 
distribution of concentration parameters of halos a t fixed mass 
is known to be a function of halo mass (e.g., INavarro et al.l 
[19971 lBullocketd]l200l . 

The baryonic component in our DMGJJR simulation follows 
the distribution of dark matter at large scales (« r~ 3 ); however, 
at small scales the gas density profile tends toward a constant 
density core that is well described by the profile of iBurkertl 
(fl995h . 
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We compute the spherically-averaged density profile for 
each of the halos in our three simulations by binning the cu- 
mulative mass profile in 100 spherical bins logarithmically 
spaced in radius from 5 - 2000/z -1 kpc centered on the dark 
matter particle with the highest local density. Dark matter 
and stellar particles are treated as point masses, while gas is 
assigned to bins by numerically integrating the gas density 
within radial annuli using 10 5 Monte Carlo samples per bin. 

We obtain fitted concentrations for each halo by minimiz- 
ing the x 1 statistic with bins weighted by the square-root of 
the bin density. In the case of particles (dark matter and stel- 
lar mass) this procedure is equivalent to assuming the bin error 
follows a Poisson distribution. Profiles are fit from an inner- 
most radius R m i n = 0.057? v j r to R w \ r . In practice the radii used 
are approximate due to the fixed widths of the bins in physi- 
cal radius. Fits to the halos in the DM0 and DMGJJR simulations 
are insensitive to the choice of innermost radius, but the den- 
sity profiles in halos in the DMG_SF simulation deviate strongly 
from the NFW profile within this radius, leading to concentra- 
tion fits that are biased to larger values as R m i n is decreased. 
The deviation is caused by the formation of a large, central 
galaxy in the DMG_SF case. In addition to fitting the total mat- 
ter distributions and the dark matter distributions for all three 
simulations, we fit the gas distributions inside the halos of the 
non-radiative simulation with the Burkert profile [Eq. ©] to 
obtain a Burkert concentration parameter c B . 

Figure [6] shows the resulting concentration parameters as 
a function of halo mass derived from these density profile 
fits. Several important features are apparent. First, consider 
the dark matter concentrations in the center panel. The ha- 
los in the DMGJJR simulation are slightly more concentrated 
(~ 10%) than the halos in the DM0 simulation, which is con- 
sistent with the effects found by p revious simulation studies 
(Rasia et al. 2004; Lin et al. 2006). As discussed earlier, dur- 
ing halo formation and relaxation the gas exchanges energy 
with the dark matter. The net result is that the gas gains energy 



while the dark matter relinquishes energy to the gas, resulting 
in slightly more compact dark matter halos and extended gas 
distributions. 

Furthermore, Figure [6] shows that the halos of the DMG_SF 
simulation are considerably more concentrated than the ha- 
los of both the DM0 and DMG_NR runs. This difference is due 
to the cooling-induced condensation of gas and formation of 
galaxies at the centers of these halos. The left panel of Fig- 
ure [6] shows the NFW concentrations that result from fitting 
the total matter distributions of the halos in each simulation. 
The results are qualitatively similar to those of the dark matter 
distributions alone, but exhibit two additional features for ha- 
los with masses <, 10 13 h~~ x M . The first feature is a relative 
decline in total matter concentration in the DMGJJR simulation 
relative to the DM0 simulation. This decline is due to the ex- 
tended, less-concentrated gas distri bution in these lower-mass 
halos, and is again consistent with lLin et al.l (120061) who find 
a smaller (~ 3%) increase in the concentration of the total 
matter distribution than for the dark matter itself. The second 
feature for low-mass halos is a dramatic rise in total matter 
concentration in the DMG_SF simulation. This increase in con- 
centration is because halos with masses < 10 13 hr l M Q cool 
gas and form galaxies more efficiently than their cluster-sized 
counterparts. The larger fraction of cold gas and stars at the 
centers of these halos gives rise to a large relative increase 
in density toward the halo centers, driving the concentration 
values higher. Finally, the right panel of Figure [6] gives the 
Burkert concentrations of the gas profiles in the DMG_NR sim- 
ulation. Note the gradual increase in c B with halo mass. This 
result will be discussed further in § 14.11 

4. INTERPRETING POWER SPECTRA WITH THE HALO MODEL 

4. 1 . Halo Model Background and Methods 

To interpret the relative differences seen in the simulated 
power spectra (Figure [2j, we turn to the halo model — an 
analytic, phenomenological framework for describing the 
clustering of dark matter, galaxies, or any other population 
associated with dark matter halos. Various aspects of the 
halo model are developed and discussed in a variety of 



studie s over several decade s dNevman & Scottlll95 2; Peebles 
1 1974 1 : iMcClelland & Sffld [l977t IScherrer & Bertschineei 
199 1[ iMa&Frvl 120001: Iseliald 12000b IScoccimarro et al 
200U IShethetalJ 12001 alibi lBerlind& Weinberg! 120021: 
Coorav & Sheth 2002). 

The halo model is predicated on the assumption that all 
of the objects of interest lie within dark matter halos. The 
two-point clustering statistics of matter are then given by the 
the sum of two terms. The "one-halo" term is from mat- 
ter elements residing within a common halo and dominates 
clustering statistics on scales smaller than typical halo sizes. 
The second "two-halo" term is from matter elements residing 
within two distinct halos and dominates on large scales. This 
decomposition is convenient because we do not expect bary- 
onic physics to alter significantly the clustering of dark mat- 
ter halos on large scales. Baryons do have significant effects 
on the structures of individual dark matter halos, as shown in 
§ 13.21 so we expect differences to be largely confined to the 
one-halo term. 

To model the total matter clustering with multiple compo- 
nents (e.g., dark matter, gas, etc.), it is conv enient to treat 
each c omponent separately, as was done by IZhan & Knoxl 
(2004). In what follows, we quote the halo model relations 
for the power spectrum. The total power is the sum of the 
power spectra of the individual components and cross terms, 



RUDD, ZENTNER, & KRAVTSOV 



20 



x 1 \ M| ^ 

\ • TOTAL 



20 





20 



L5 



10 1; 



10 13 
M [h- 



10 1 ' 



I 



Mo 



- 


GAS 


- 

: y/// ; i - 


//>.. 

A// Ay 
/ /I 






1 III 


M 1,1 


12 10 13 


10" 


M [h-> M Q ] 





^(&) = Ei,j fifjPijik), where /; = f2,/Q M refers to the univer- 
sal mass fraction in the ith matter component. The one-halo 
contribution is given by 



1 

Pipj 



2 dn 
Am m fj(m)fj(m) — — Ai(k; m)Aj(k; to), 



(5) 



where p, is the mean density in the ith matter component, 
fi(m) is the average fraction of mass in halos of total mass 
to residing in the ith component, An I dm is the mass function 
of halos, and A,(k; to) is the Fourier transform of the mean 
density profile of the ith component in halos of total mass 
to. For example, the profiles of dark matter halos are of- 
ten modeled by NFW profiles, in which case, Ai(k; to) is the 
Fourier transform of the NFW density profile (e.g., given by 
IScoccimarro et al J 120011) with a concentratio n parameter set 
by some relation (e.g., Bull ock et al.ll200ll) . The two-halo 
contribution to P(k) is 



Fig. 6. — Mean concentration of radial distributions of main mass components as a function of halo mass. Shaded regions contain 68% of halos in a given mass 
bin. Scatter in the concentration relation for the DM0 simulation is comparable to that in the DMG.NR simulation and is omitted for clarity. Error bars correspond 
to the estimated error on the mean concentration in each bin. Left panel: the mean NFW concentration fit to the total (i.e., DM+gas+stars) mass distribution for 
all three of our simulations as indicated at the bottom of the panel. Center panel: the mean NFW concentration fit to the dark matter distribution in all three 
simulations. Right panel: the mean Burkert concentration fit to the gas density profiles in the DMG_NR simulation. 

enforced. As halos have a finite extent set by their virial 
radii, the integrals in Eq. (0 should not extend over all mass 
but should be limited to halos with virial radii smaller than 
r ~ k~ x . This effect is known as halo exclusion. Though more 
complex and accurate imp lementations of halo exclusion ex- 
ist (e.g., iTinker et al.ll2.Q06h, we u se the model for halo exclu- 
sion introduced by IZheng d2004l) . Briefly, we set the upper 
bounds on the integrals in Eq. (O to the halo mass that cor- 
responds to a virial radius of r max = 2nlC l . Previous studies 
have found this prescrip t ion to be useful for practical applica- 
tions (e.g., [Zhin3200l|Zehayi^ 

The last ingredients necessary to build a halo model of the 
matter power spectrum are specifications of the density pro- 
files that characterize the distribution of each matter compo- 
nent within halos. We treat each of the cases of pure dark 
matter, dark matter with non-radiative gas, and dark matter 
with gas cooling and star formation slightly differently, with 
prescriptions motivated by our set of simulations. 

We model the dark matter halos in both the A-body and 
non-radiative cases with the NFW density profile [Eq. (fj)]. 
As in § 13.2.21 the concentrations of halos are different in each 
case, and we include this effect in our implementation of the 
halo model. In our modeling, it is necessary to extrapolate be- 
yond the range of concentrations probed directly by our sim- 
ulations. Partly motivated by the fact that we aim to represent 
the features of our simulated spectra qualitatively, we adopt a 
particular for m of the analy t ic mod el for halo concentrations 
introd uced in | Bullock et alj d2001[). Similar t o other authors 
(e.g. iDolag et all l2004t iKuhlen et al J 120051: [Wechsler et alj 
12006b iMaccio etal.1 120061) . we find that the relationship be- 
tween concentration and mass in our simulations has a smaller 
normalization and a slightly shallower slope than that of the 
iBullock et al.l d2001l) model in its original form. We find that 
the mean concentration as a fun ction of mass in the D M0 sim- 
ulation is well described by the Bullo cket alj ((2001) model 
with parameters F = 10~ 5 and K = 1.7. We stress that these 
parameters are not the result of a formal fitting procedure and 
defer further exploration of the concentration-mass relation to 
future work. 

The halos in our DMGJJR simulation exhibit somewhat 
higher concentrations than those in the DM0 simulation. Over 
the mass range measured in the simulation, we use the mea- 
sured concentration-mass relation (e.g., Fig. [6]) from the sim- 
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P lm (k) is the linear matter power spectrum, and bh(m) is the 
mass-dependent halo bias. 

Our primary aim in applying the halo model is to study the 
qualitative features of the spectra from our simulations rather 
than to provide a precise, quantitative description. There- 
fore, we adopt the fitting forms for the m ass function and 
linear bias of dark matter halos provided by Sheth & Tormen 
(1999), rath er than any of several updated bias pre scrip- 
tions (e.g., Uenkins et alj 1200 U ISeliak & Warrenl [2004 see 
ICoorav & Shethl2002l and lz"entnerl2006l recent reviews). This 
choice guarantees that the two normalization relations 
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are satisfied identically without making any further, and of- 
ten arbitrary, choices about how these relations should be 
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ulations and use the Bull ock et al.l d2001l) -like model in order 
to extrapolate outside of the ran ge of measured halo concen- 
trations. Specifically, we fix the lBullock etai] (1200 ll) param- 
eter F = 10~ 5 as for the dissipationless simulation and al- 
low K to float in order to best match the normalization of the 
concentration-redshift relation at the redshift of each simula- 
tion output. Typical values of K are ~ 10% larger than in the 
dissipationless case. 

Including gas in the halo model for the non-radiative simu- 
lation introduces additional complexity into the standard halo 
model. It is now necessary to specify how much gas there 
is inside each halo. Figure [5] gives the mean fraction of halo 
mass in baryons as a function of halo mass for both simu- 
lations with baryons. In the non-radiative case we adopt a 
constant value of /b = 0.94£2b/£2m for all halos. This value 
is consistent with the baryon fractions in the largest halos of 
DMGJJR to redshifts greater than z — 3. Treating the baryon 
fraction in this way is sensible because the one-halo contribu- 
tion to the 3D and convergence power spectra should be dom- 
inated by the largest halos on scales of interest (see Fig. 4 of 
IZhan & Knoxl (|2004)) and insures that we do not model a res- 
olution effect. For simplicity, we treat the remaining baryons 
as an unclustered component. 

Our non-radiative simulation revealed that the NFW pro- 
file is not a good description of the distribution of gas within 
halos (§ I3.2.21 i. In this case we model the gas separately 
with a Burkert profile. The constant-density core within 
r <k /?vir/cB causes a suppression in the small-scale power 
at high wavenumbers relative to an NFW profile. The trend of 
cb as a function of halo mass defined in this way is the oppo- 
site of the NFW concentrations of halos. An exploration of the 
possible source of this trend is beyond the scope of this work. 
We interpolate the mean concentration relation in Figure [6] 
linearly in log(cg) as a function of log(m). For masses below 
our resolution limit we assume a constant eg equal to that of 
the lowest mass bin, although we have explicitly checked that 
our results are insensitive to this choice. That our choice of 
extrapolation does not affect our results significantly is unsur- 
prising, as halos below 10 12 hr l M contribute only a small 
fraction of the power on scales of interest. Likewise, we set 
cb equal to the value at the highest mass bin for halos larger 
than the halos in our simulation volume. 

In principle, modeling power spectra in the case where ra- 
diative cooling and star formation are included is more com- 
plex still. The halos in our simulations with radiative cooling 
and star formation exhibit baryon density profiles that are sig- 
nificantly more complex than those in the non-radiative simu- 
lations. Instead of a single, hot baryonic component, cooling 
and star formation give rise to a condensed component of cold 
gas at the center of each halo, a fraction of which turns into 
stars. Moreover, the baryons in the gas phase separate into 
hot and cold components that do not follow a common den- 
sity profile. 

Describing all three components of the baryons separately 
would be difficult and, while interesting, we relegate this 
to future work focusing on galaxy formation physics. For- 
tunately, such a detailed exploration of the distribution of 
baryons within halos is not necessary in order to understand 
the qualitative features in the relative power spectra. The pro- 
cess that causes the strong condensation of baryons at halo 
centers also affects the dark matter, and the total matter profile 
is well described by an NFW profile outside as 0.05 - O.LRvir- 
As such, we model the total mass profiles of halos in the case 
of radiative cooling and star formation as a single compo- 
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Fig. 7. — The halo model power spectra for the non-radiative (DMG_NR) and 
cooling+star formation (DMG.SF) cases are plotted in left and right columns 
respectively. Line types and panel layout follow the convention of Figure [2] 
The dot-dashed lines represent the one- and two-halo contributions to the 
baryon power spectra (for the DMG.SF case baryons are assumed to trace the 
total matter distribution). The two-halo term dominates at low k, while the 
one-halo term is considerable at intermediate wavenumbers. Note that the 
simple modifications to the standard halo model formalism described in the 
text reproduce the qualitative features of the simulated spectra well. 

nent described by an NFW profile with a concentration set 
according to the tabulated mean concentration-mass relation 
extracted from the cooling simulations (left panel of Fig. [6j. 
We extrapolate outside of the mass range probed by our sim- 
ulations in the same manner as for dark matter halos in the 
non-radiative case. 

4.2. Halo Model Power Spectra 

The left-column of Figure [7] shows the power spectra from 
a halo model including the effects of non-radiative baryonic 
physics as described above. At z — there is good qualita- 
tive agreement between the halo model power spectra and the 
power spectra from the simulations. In particular, the qual- 
itative features of an enhancement in the dark matter spec- 
trum at high wavenumber (k > 1/iMpc -1 ), a comparably sup- 
pressed gas power spectrum, and a total matter power spec- 
trum with an enhancement at intermediate wavenumber (k ~ 
1 - lOh Mpc -1 ) and subsequent decrease at high wavenumber 
(k > 20/iMpc -1 ) are all reproduced by this simple, augmented 
halo model. 

All of these features are set by modifications to the one-halo 
term in the halo model. The increase in the dark matter power 
is caused by the increased concentrations of halos in the non- 
radiative case. The relative decline of the gas power spectrum 
is due to the more extended profile of the gas relative to an 
NFW profile [see Eq. (@)]. The two features in the total mass 
power spectrum are a reflection of both of these effects. 

At higher redshifts these features move to larger wavenum- 
bers, a shift that is also present in the simulations but to a 
lesser degree. In the context of the halo model, this shift 
in scale is caused by the change of the peak location of the 
one-halo term as the sizes of the largest halos decrease with 
increasing redshift. We attribute differences between the red- 
shift evolution of the simulations and of the halo model to 
differences in the evolution of the halo mass function and bias 
relations in the simulation, which are not exactly given by 
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the analytic fSheth & Tormenl ( 1 19991) relations. In particular, 
the mass function exhibits large variations at the high-mass 
end due to the limited simulation volume. There also exists 
an apparent large-scale bias of the gas component because 
some fraction of the gas in our simulations is not contained 
within halos and we have treated this fraction as an unclus- 
tered component, whereas in the simulations these baryons 
are clustered, albeit somewhat less strongly than the baryons 
inside halos. 

The right-column of Figure [7] shows the halo power spec- 
tra resulting from our model of the radiative case. Again, the 
qualitative features of the simulated spectra are reproduced by 
the halo model. The power increases steeply at k > 1/zMpc ~' 
at all redshifts, although the transition is quicker at lower 
redshift, reaching a ~ 40% increase of total power by k ~ 
4-h Mpc 1 at z — 0. At higher redshifts this transition shifts to 
smaller scales. Again, the relative importance of the one-halo 
term is responsible for setting the scale that this occurs. 

The increase of power in the one-halo term is caused pri- 
marily by the larger halo concentrations in the DMG_SF run. 
Increasing the normalization of the c(m) relation to the level 
observed in this simulation leads to a considerable enhance- 
ment of the matter power spectrum at small scales. Finally, 
note that the slopes of the simulated power spectra continue 
to increase as a function of wavenumber at small scales, a fea- 
ture not present in our halo model spectra. This trend is likely 
due to the fact that our simple halo model is not designed to 
reflect the contributions of the condensed gas and stellar com- 
ponents within r < 0.05 - 0.17? v j, of each halo, deviations 
from the NFW form, or halo substructure, and so it cannot 
represent the matter power spectra at very high wavenumbers 
well. 

4.3. An Heuristic Model For Power Spectra With Cooling 
and Star Formation 

The trends discussed in the previous subsections suggest 
that the primary effect of baryons is to redistribute matter 
toward the centers of halos. In the case where the baryon 
component is allowed to cool, this effect is particularly dra- 
matic, and can be adequately described over a wide range of 
scales by an increase in the effective NFW concentration (see 
§ |3.2.2t . Our results indicate that the differences are largely 
relegated to the one-halo term in the halo model [Eq. ©]. 

This physical intuition suggests a simple method to model 
the effect of baryonic cooling on the power spectrum over 
the limited range of wavenumbers important for weak lensing 
studies. The Fourier transform of the NFW profile is analytic 
and given by (e.g., [Scoccimarro et al. 2 001b 



1 



S in(77)[Si([l+c]77)-Si(77)] 
Ci(77)] 



/nfw(c) 

+ cos(77)[Ci([l + c];/) 
sin(?7) 
[1 + c]t) 



(10) 



where c is the halo concentration, 77 = kR v i r /c, R vu is the halo 
virial radius, /nfwM = ln(l + x) - x/(l + x), and Si(x) and 
Ci(x) are the sine and cosine integrals, respectively. Eq. ( TTOb 
is normalized such that the integral of A(x) over all space is 
unity. The systematic difference in halo profiles as a function 
of mass is encoded in the shift in halo concentrations at that 
mass. The power spectrum on scales of order k ~ Ih Mpc 1 
is dominated by the one-halo term contribution from cluster- 



mass halos (~ 10 h M Q ). As such, we model the modi- 
fication to the power spectrum due to cooling on scales near 
a few h Mpc 1 by multiplying the dissipationless power spec- 
trum by a ratio of Fourier-transformed NFW profiles, 



Ps F (k) * PdmW 



A(R vil -k/c 2 ,C2) 



A{R vir k/c u ci) 



(11) 



Roughly speaking, the second factor on the right side of 
Eq. ( fTTT l represents a shift in the power spectrum due to 
a boost in the concentrations of halos which dominate the 
largest nonlinear scales. This factor is unity on scales larger 
than the typical scales of halos k < R'l and differs from unity 
at higher wavenumbers as one would expect for a one-halo 
effect. To implement this model, one chooses a characteris- 
tic R vlr , a concentration c\ that represents the halos in a typi- 
cal dark matter simulation, and a second concentration q that 
represents the effective concentrations of the mass profiles of 
halos with baryonic cooling and galaxy formation included. 

While halo properties evolve, in our model of the conver- 
gence power spectrum the lensing weight function has a broad 
peak centered near z w 0.55, due to our choice of a single thin 
sheet of sources at z s = 1. We therefore select parameters for 
our model to represent typical halos at this redshift, neglecting 
possible evolution. At z = 0.55, the largest halos in our sim- 
ulations have masses ~ 10 14 hr x M Q , corresponding to virial 
radii of R v i r ~ Lift -1 Mpc comoving. In the dissipationless 
simulation, DM0, halos of this mass have mean concentrations 
C\ ~ 5. The concentrations in the DMG_SF simulation are en- 
hanced by a factor of a 1.7 and so we take c 2 - 1.7ci =8.5. 

Figure [8] shows the three-dimensional and convergence 
power spectra that result from applying our model with these 
parameters to the power spectra obtained from the DM0 simu- 
lation. Notice several aspects of Figure [8] which have changed 
relative to Figure [2] and Figure [3] First, we have removed the 
large-scale offset between the power spectra of the DMG_SF 
and DM0 simulations by forcing them to be equivalent at the 
lowest wavenumber, k w 0.13/; Mpc -1 . Second, in the right 
panel of Figure [8] we now plot convergence power spectra 
relative to the spectrum of the DMG_SF simulation rather than 
relative to the DM0 spectrum as in Fig. [3] The DM0 power 
spectrum now runs downward toward negative residual in- 
stead of the residual of the cooling run trending upward with 
wavenumber. For comparison, the dotted line in the left panel 
of Figure [8] shows the result of adding a Poisson term to the 
dissipationless power spectrum as an alternative model for the 
net effect of cooling. The motivation for a Poisson model is to 
treat the compact halo cores as point masses. Note that the 
Poisson model does not reproduce the correct power spec- 
trum shape, regardless of amplitude. Moreover, the large 
number densities required to produce even reasonable agree- 
ment (~ 0.1 Z; 3 Mpc -3 ) in this regime correspond to halos of 
mass ~ 10 11 hT x M Q which is significantly smaller than the 
halos that dominate the power on the range of scales we con- 
sider, and thus the original physical motivation for the Poisson 
model is not substantiated. 

The model of enhanced NFW concentrations, on the other 
hand, successfully tracks the departure of the DMG_SF power 
spectrum at scales I < 4000 from the DM0 spectrum. Thus 
we extend the range of I where we can accurately model 
the power spectrum by an order of magnitude at a cost of 
three extra parameters. The model fails on small scales 
I > 10 4 for at least two reasons. The first is that the distri- 
butions of cold gas and stellar mass are not well described 
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Fig. 8. — A model for the gross effect of baryon cooling on the matter power spectrum. The left panel shows matter power spectra at z = 0.55 (where the 
lensing weight function peaks given a thin plane of sources at redshift z s = 1), relative to the matter power spectrum in the A'-body simulation. In this case, we 
have scaled out the small, large-scale bias of the power in the cooling simulation relative to the A'-body simulation at small wavenumber. The solid line shows 
the result of the cooling simulation. The dotted line shows an example of adding a Poisson term to the power spectrum of the A'-body simulation. The dashed 
line shows the result of taking our simple model for the net effect of baryonic dissipation on the matter power spectrum. In this particular example, we take 
i? v j r = \.\h Mpc, C[ = 5, and ci = \.lc\ = 8.5. The right panel shows convergence power spectra relative to the power spectrum in the DMG.SF simulation. 
Note that this choice of normalization is different from FigurefS] which is why why the DM0 case here appears as a deficit. The line types are the same as in the 
left panel and the error bands are the same as those in Figure[5] 



by the high-concentration NFW profile prescription within 
ss 0.05 - 0.1/? v i r , so the model should break down when the 
central regions of halos contribute significantly. The second 
is that our choice of parameters for this modification is driven 
by the properties of the largest halos in our simulations. This 
choice of parameters should not continue to produce a viable 
model on scales where the structure of smaller halos or of 
subhalos contributes significantly to the power. 

Our fiducial choice of parameters R w \ r , c\, and q reflects 
the properties of the largest halos in the simulation volume 
at redshift z = 0.55. In practice one should probably intro- 
duce these parameters or a subset thereof (for example, one 
could suppose that the concentrations of halos in dissipation- 
less simulations are well understood, fixing ci) in any data 
analysis and marginalize over our uncertainty in these nui- 
sance parameters. Ideally one would account for the evolution 
of these parameters with redshift, perhaps by parameterizing 
this evolution or by introducing independent sets of parame- 
ters for different tomographic redshift bins. It remains to be 
determined how much of the constraining power from high-f 
modes will be lost in this marginalization process. We expli- 
cate the evolution of these parameters and the degradation of 
cosmological parameters arising from the additional nuisance 
parameters in a forthcoming study. 

Our results suggest that modifications to the power spec- 
trum are dominated by the one-halo term, which motivates 
this simple model. We have not attempted a more quantita- 
tive analysis because it is unclear that such an analysis is war- 
ranted given the intrinsic limitations of our simulation data. 
Clearly, further theoretical efforts are needed to constrain and 
understand the effects of galaxy formation physics and refine 
the analytic model of these effects. A first possible refine- 
ment would be to scale dissipationless spectra by the ratio 
of two complete one-halo terms [the form of which is given 
in Eq. ©]. In the upper term, the halos could follow a dif- 
ferent concentration-mass relation than dissipationless halos 
(for example, parameterized by a different power law) in or- 



der to model the modified structures of halos with baryonic 
effects included. Another possible refinement (at the cost of 
additional nuisance parameters) would be to add additional 
components to the power spectrum modeling to account for 
the distributions of hot and cold gas and stars within halos 
and subhalos or to allow for more general profile shapes than 
NFW. 

5. DISCUSSION AND CONCLUSIONS 

We have performed a theoretical study of the impact of 
baryonic physics on predictions of the matter power spectrum 
as it will be measured by contemporary and forthcoming weak 
lensing surveys. With the exception of the recent study by 
iJing et all (2006), such predictions have hitherto been derived 
from dissipationless Af-body simulations that neglect the ad- 
ditional physics of the baryonic component of the universe, 
or by semi-analytic models (e.g., IZhan & Kn ox 2004; IWhitd 
2004) which do not treat baryons self-consistently. Yet it is 
of paramount import ance that these predicti ons be extremely 
precise (< 1%, e.g. lHuterer & Taka da 2005) in order to real- 
ize fully the power of weak lensing surveys to constrain the 
properties of dark energy. 

We studied these effects using a suite of numerical simu- 
lations of cosmological structure formation. Each simulation 
started from the same initial density field but differed in the 
treatment of baryons in the universe. We also used the ana- 
lytic and phenomenological halo model of clustering in order 
to interpret our results and provide guidance for future stud- 
ies. Our primary results can be summarized as follows. 

The effect of baryons on the matter power spectrum is sig- 
nificant at a level that may be important for the interpretation 
of contemporary weak lensing surveys and will almost cer- 
tainly need to be accounted for in the interpretation of future 
surveys. The enhancement in power in the DMG_SF simulation 
relative to the DM0 simulation due to baryon dissipation and 
star formation at £ > 800 is large enough that it may affect the 
interpretation of the results from contemporary surveys like 
the CFHT Wide Survey. This enhancement is much larger 
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than the statistical error of future LSST- and SNAP-like sur- 
veys on scales that have been used in previous forecasts of the 
constraining power of such surveys {{ < 3 x 10 3 ). We recog- 
nize that the stellar masses of the galaxies in our DMG_SF sim- 
ulation are large compared to observed galaxies and take this 
to indicate that the enhancement that we quote for the DMG_SF 
simulation is likely to be an overestimate. Nevertheless, we 
expect these effects to be quite important for future surveys. 
Even in the absence of cooling, future weak lensing surveys 
would still be sensitive to the effects of different dynamics of 
dark matter and baryons and their different final distributions 
in virialized halos. 

Our findings are in qualitative agreement with the results 
of iJing et al.l d2()()6). However, there are quantitative differ- 
ences, which we attribute to different numerical hydrodynam- 
ics methods, and different prescriptions for star formation and 
stellar feedback. These differences serve to emphasize the fact 
that power spectrum predictions are sensitive to the poorly- 
understood physics of galaxy formation and the specific de- 
tails of the treatments of the baryonic components in different 
cosmological simulations. 

The effects of baryons on the power spectrum on large 
scales (k < 1/jMpc" 1 or t < 800) are minimal. In fact, we 
find that most of the identified effects on P(k) are explained 
by changes in the mass distributions within virialized objects 
in response to baryonic physics. In the simulation which in- 
cludes radiative cooling, large amounts of gas cool and con- 
dense at halo centers where they form stars. As the baryons 
cool to halo centers, they drag dark matter with them. The 
net result is a large enhancement in the relative concentration 
of matter toward the centers of halos resulting in a dramatic 
increase in small-scale power. On the other hand, the gas in 
the DMGJJR simulations is unable to dissipate energy gained 
during gravitational collapse. Dark matter halos in this case 
have their concentrations slightly enhanced leading to an in- 
crease in power on intermediate scales (1 < fc//iMpc~' < 20 
or 800 < i < 10 4 ), while the baryons associated with halos 
remain hot and extended leading to a net deficit in power on 
small scales (k/hMpc^ > 20 or t > 10 4 ). We have used 
detailed halo models of the matter power spectra designed to 
mimic our simulation results in order to support this interpre- 
tation. 

Observational measurements of the cluster concentration- 
mass relation would be of significant help in constraining 
implementations of galaxy formation physics in simulations 
and in testing the overall adequacy of dissipationless simula- 
tions. Such observational constraints can be derived from ei- 
ther weak or strong lensing, or through X-ray me asurements 
of the intra-cluster medium (see f or example, Buo te et all 
l2006t IComerford & Nataraianll2007h . At present, these mea- 
surements have their own sets of assumptions and system- 
atic uncertainties and their interpretation requires further tests 
against numerical simulations. However, we envision that 
such constraints can become significantly tighter in the future. 

One of the primary points of this paper is that it is imper- 
ative that a discussion begin in earnest in order to determine 
how best to deal with the important, yet poorly-understood 
physics of galaxy formation in the context of forthcoming 
weak lensing surveys. In § 14.31 we showed how the halo 
model of clustering could be used to model the basic features 
of the matter power spectrum in the presence of baryons quite 
successfully and at a relatively minor cost. The model relies 
on a parameterized form for the one-halo term in P(k), the 



parameters of which would be fit to either simulation or ob- 
servational data as it becomes available. We stress that the 
modifications proposed in § 14.31 work well to describe our 
simulations data, but we have not yet tested this approach at 
the requisite level of detail and on a wide enough spectrum 
of simulation results to have confidence that it will describe 
observational data precisely. Further, we have not studied to 
what degree the additional parameters required by this effec- 
tive, parametric model of the effect of baryons on the matter 
power spectrum will degrade future constraints on dark en- 
ergy parameters. We plan to attack both of these issues in a 
forthcoming study. 

As for the prospects for direct predictions of the matter 
power spectrum at k > l/iMpc~' or I > 10 3 , the chal- 
lenge facing theorists is and will remain significant. Simu- 
lations which treat the effects of baryons and galaxy forma- 
tion physics are currently limited in dynamic range. As a 
consequence, simulations that resolve small scales and fol- 
low processes relevant for galaxy formation are unable to 
treat volumes that are large enough to render cosmic vari- 
ance insignificant and even these state-of-the-art simulations 
make assumptions about the net effects of processes that oc- 
cur on scales that are unresolved. Unlike dissipationless N- 
body simulations, the computational cost of simulations that 
include these additional physical processes is large enough 
that a direct exploration of the way in which baryonic effects 
change with cosmology will be impractical for some time. 

Additionally, it is unclear whether current simulations with 
galaxy formation have reached numerical convergence at the 
present achievable resolution and with the included physics. 
For example, current simulations produce stellar fractions in 
galaxy cluster halos that are a factor of two to three too large 
compared with observational estimates of » 15% — 20% of 
the un iversal baryon fraction dLin et al.ll2.Q03l: iGonzalez et al.l 
120071) . Some of the tension may be alleviated if a signifi- 
cant fraction of stars are in a diffuse co mponent not normally 
accounted for in observa tions (e.g., see lGonzalez et al.ll2005l 
l2007t ISeigar et al1l2006l) . There is additional, indirect evi- 
dence from observations of the intracluster gas via X-rays and 
Sunyaev-Zel'dovich effects that the amount of baryons con- 
densi ng out of the hot gas is approximately co rrect in simula- 
tions (lAfshordi et al.ll2006t iNagai et al.ll2007l) . Still, the role 
of feedback from Active Galactic Nuclei or additional plasma 
physics in halo cores and their effects on the distribution of 
baryons and the power spectrum are also yet to be understood. 

In summary, ab initio predictions for the matter power spec- 
trum as a function of cosmology on relevant scales are not 
feasible, so alternative strategies must be explored. We sug- 
gest that a worthwhile, general approach is to attempt to iso- 
late key effects of galaxy formation and encapsulate them in 
as few parameters as possible. Such parameterized analytic 
models could be used in conjunction with predictions of large- 
volume, dissipationless simulations to model the observable 
convergence power spectra at the expense of a small set of 
nuisance parameters. As we have shown in this work, this 
approach appears to be quite promising if the primary effect 
of galaxy formation processes relevant to convergence power 
spectra is to modify the distribution of mass within halos. 

We would like to thank Hu Zhan for useful discussions on 
the possible effects of baryon dissipation on the power spec- 
trum during the early stage of this study and for helpful com- 
ments on an early version of this manuscript. We are grate- 



BARYONS AND THE POWER SPECTRUM 



13 



ful to Wayne Hu, Dragan Huterer, Jeremy Tinker, and Mar- 
tin White for many useful discussions throughout the course 
of this project, to Sarah Hansen for several careful readings 
of the manuscript, and to Anatoly Klypin for generating and 
providing initial conditions and snapshots of the dissipation- 
less simulation. This project was supported by the National 
Science Foundation (NSF) under grants No. AST-0239759 
and AST-0507666, by NASA through grant NAG5- 13274, 
and by the Kavli Institute for Cosmological Physics (KICP) 
at the University of Chicago. ARZ is supported by The NSF 



Afshordi, N., Lin, Y.-T., Nagai, P., & Sanderson, A. J. R. 2006. 

MNRAS submitted I astro-ph/06 1270(5} 
Albrecht, A., Bernstein, G., Cahn, R., Freedman, W. L., Hewitt, J., Hu, W., 

Huth, J., Kamionkowski, M., Kolb, E. W., Knox, L., Mather, J. C, Staggs, 

S., & Suntzeff, N. B. 2006, 1 astro-ph/0609591) 
Aldering, G. 2005, New Astronomy Review, 49, 346 

Aldering, G., Althouse, W., Amanullah, R., Annis, J., & et al. 2004, PASP 

submitted ( astro-ph/0405232 1 
Annis, J., Castander, F. J., Evrard, A. E., & et al. 2005, lastro-ph/0510194 1 
Astier, P., Guy, J., Regnault, N., & et al. 2006, A&A, 447731 
Balogh, M. L., Pearce, F. R., Bower, R. G., & Kay, S. T. 2001, MNRAS, 326, 

1228 

Berlind, A. A. & Weinberg, D. H. 2002, ApJ, 575, 587 
Blumenthal, G. R., Faber, S. M., Flores, R., & Primack, J. R. 1986, ApJ, 301, 
27 

Borgani, S., Governato, E, Wadsley, J., Menci, N, Tozzi, P., Quinn, T., Stadel, 

J., & Lake, G. 2002, MNRAS, 336, 409 
Bryan, G. L. & Norman, M. L. 1998, ApJ, 495, 80 

Bullock, J. S., Kolatt, T. S., Sigad, Y, Somerville, R. S., Kravtsov, A. V., 
Klypin, A. A., Primack, J. R., & Dekel, A. 2001, MNRAS, 321, 559 

Buote, D. A., Gastaldello, E, Humphrey, P. J., Zappacosta, L., Bullock, J. S., 
Brighenti, E, & Mathews, W. G. 2006, lastro-ph/0610135 I 

Burkert, A. 1995, ApJ, 447, L25+ 

Comerford, J. M. & Natarajan, P. 2007, 1 astro-ph/0703 126 1 
Cooray, A. & Hu, W. 2001, ApJ, 554, 56 
Cooray, A. & Sheth, R. 2002, Phys. Rep., 372, 1 

Ci-ain, R. A., Eke, V. R., Fren k, C. S„ Jenkins, A ., McCarthy, I. G., Navarro, 

J. E, & Pearce, F. R. 2006, 1 astro-ph/06 10602 1 
Dave, R., Cen, R., Ostriker, J. P., Bryan, G. L., Hernquist, L., Katz, N., 

Weinberg, D. H., Norman, M. L., & O'Shea, B. 2001, ApJ, 552, 473 
David, L. P., Jones, C, & Forman, W. 1995, ApJ, 445, 578 
Dodelson, S. & Zhang, P. 2005, Phys. Rev. D, 72, 083001 
Dolag, K, Bartelmann, M., Perrotta, E, Baccigalupi, C, Moscardini, L., 

Meneghetti, M., & Tormen, G. 2004, A&A 
Eisenstein, D. J. & Hu, W. 1999, ApJ, 511, 5 

Eisenstein, D. J., Zehavi, I., Hogg, D. W., & et al. 2005, ApJ, 633, 560 
Eke, V. R., Navarro, J. E, & Frenk, C. S. 1998, ApJ, 503, 569 
Ettori, S., Dolag, K., Borgani, S., & Murante, G. 2006, MNRAS, 365, 1021 
Ferland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. W., Kingdon, J. B., 

& Verner, E. M. 1998, PASP, 1 10, 761 
Frenk, C. S., White, S. D. M., Bode, P., & et al. 1999, ApJ, 525, 554 
Gnedin, O. Y, Kravtsov, A. V., Klypin, A. A., & Nagai, D. 2004, ApJ, 616, 

16 

Gonzalez, A., Zaritsky, D., & Zabludoff, A. 2007, ApJ submitted 

Gonzalez, A. H., Zabludoff, A. I., & Zaritsky, D. 2005, A pJ, 618, 195 

Gottloeber, S., Yepes, G., Wagner, C., & Sevilla, R. 2006, 1 astro-ph/0608289 1 

Haardt, F. & Madau, P. 1996, ApJ, 461, 20 

Hagan, B., Ma, C.-P, & Kravtsov, A. V. 2005, ApJ, 633, 537 

Heavens, A. 2003, MNRAS, 343, 1327 

Heitmann, K, Ricker, P. M., Warren, M. S., & Habib, S. 2005, ApJS, 160, 28 
Hockney, R. W. & Eastwood, J. W. 1981, Computer Simulation Using 

Particles (Computer Simulation Using Particles, New York: McGraw-Hill, 

1981) 

Hu, W. 1999, ApJ, 522, L21 

Hu, W. & Tegmark, M. 1999, ApJ, 514, L65 

Huterer, D. 2002, Phys. Rev. D, 65, 063001 

Huterer, D. & Takada, M. 2005, Astroparticle Physics, 23, 369 

Huterer, D., Takada, M., Bernstein, G., & Jain, B. 2006, MNRAS, 366, 101 

Jenkins, A., Frenk, C. S., Pearce, F. R., Thomas, P. A., Colberg, J. M., White, 

S. D. M., Couchman, H. M. P., Peacock, J. A., Efstathiou, G., & Nelson, 

A. H. 1998, ApJ, 499, 20 
Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., Cole, S., Evrard, 

A. E., Couchman, H. M. P., & Yoshida, N. 2001, MNRAS, 321, 372 



through the Astronomy and Astrophysics Postdoctoral Fel- 
lowship Program under grant AST-0602 1 22 and by The KICP 
at The University of Chicago. The cosmological simulations 
used in this study were performed on the SGI Altix sys- 
tem (Columbia) at NASA Ames and the IBM RS/6000 SP4 
system (copper) at the National Center for Supercomputing 
Applications (NCSA). In addition, we have made extensive 
use of the NASA Astrophysics Data System and arXiv.org 
preprint server. 



Jing, Y. P., Zhang, P., Lin, W. P., Gao, L., & Springel, V. 2006. ApJ, 640, 
L119 

Kaiser, N. 1992, ApJ, 388, 272 
— . 1998, ApJ, 498, 26 

Katz, N. & White, S. D. M. 1993, ApJ, 412, 455 
Kennicutt, R. C. 1998, ApJ, 498, 541 

Klypin, A., Gottlober, S., Kravtsov, A. V, & Khokhlov, A. M. 1999, ApJ, 
516, 530 

Kravtsov, A. V. 1999, Ph.D. Thesis 

Kravtsov, A. V., Klypin, A., & Hoffman, Y. 2002, ApJ, 571, 563 
Kravtsov, A. V. & Klypin, A. A. 1999, ApJ, 520, 437 
Kravtsov, A. V., Klypin, A. A., & Khokhlov, A. M. 1997, ApJS, 1 1 1, 73 
Kravtsov, A. V. et al. 2004, ApJ, 609, 35 

Kuhlen, M., Strigari, L. E., Zentner, A. R., Bullock, J. S., & Primack, J. R. 

2005, MNRAS, 357, 387 
Levine, R. & Gnedin, N. Y. 2006, ApJ, 649, L57 

Lewis, G. E, Babul, A., Katz, N., Quinn, T, Hernquist, L., & Weinberg, D. H. 

2000, ApJ, 536, 623 
Limber, D. N. 1954, ApJ, 119, 655 

Li n, W. P., J ing, Y. P., Mao, S., Gao, L., & McCarthy, I. G. 2006, 

I astro-ph/0607555 I 
Lin, Y.-T., Mohr, J. j', & Stanford, S. A. 2003, ApJ, 591, 749 
Ma, C.-P. & Fry, J. N. 2000, ApJ, 543, 503 
Ma, Z., Hu, W., & Huterer, D. 2006, ApJ, 636, 21 

Maccio, A. V., Dutton, A. A., van den Bosch, F. C, Moore, B., Potter, D., & 

Stadel, J. 2006, MNRAS submitted I astro-ph/0608 157 I 
McClelland, J. & Silk, J. 1977, ApJ, 217, 331 
Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007, ApJ submitted 
Navarro, J. E, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Neyman, J. & Scott, E. L. 1952, ApJ, 116, 144 

Pearce, F. R., Thomas, P. A., Couchman, H. M. P., & Edge, A. C. 2000, 

MNRAS, 317, 1029 
Peebles, P. J. E. 1974, ApJ, 189, L51 

Perlmutter, S., Aldering, G., Goldhaber, G., & et al. 1999, ApJ, 517, 565 
Rasia, E., Tormen, G., & Moscardini, L. 2004, MNRAS, 351, 237 
Refregier, A. 2003, ARA&A, 41, 645 

Refregier, A., Massey, R., Rhodes, J., Ellis, R., Albert, J., Bacon, D., 
Bernstein, G., McKay, T, & Perlmutter, S. 2004, AJ, 127, 3102 

Riess, A. G., Filippenko, A. V., Challis, P., & et al. 1998, AJ, 116, 1009 

Riess, A. G., Strolger, L.-G., Tonry, J., Casertano, S., Ferguson, H. C, 
Mobasher, B., Challis, P., Filippenko, A. V., Jha, S., Li, W., Chornock, 
R., Kirshner, R. P., Leibundgut, B., Dickinson, M., Livio, M., Giavalisco, 
M., Steidel, C. C, Bem'tez, T., & Tsvetanov, Z. 2004, ApJ, 607, 665 

Scherrer, R. J. & Bertschinger, E. 1991, ApJ, 381, 349 

Scoccimarro, R., Sheth, R. K., Hui, L., & Jain, B. 2001, ApJ, 546, 20 

Se igar, M. S„ Grah am, A. W., & Jerjen, H. 2006, MNRAS submitted 
I astro-ph/06 12229 1 

Seljak, U. 2000, MNRAS, 318, 203 

Seljak, U. & Warren, M. S. 2004, MNRAS, 355, 129 

Sellwood, J. A. & McGaugh, S. S. 2005, ApJ, 634, 70 

Semboloni, E., van Waerbeke, L., Heymans, C, Hamana, T, Colombi, S., 
White, M., & Mellier, Y. 2006, MNRAS submitted I astro-ph/0606648 I 

Sheth, R. K, Diaferio, A., Hui, L., & Scoccimarro, R. 2001a, MNRAS, 326, 
463 

Sheth, R. K, Hui, L., Diaferio, A., & Scoccimarro, R. 2001b, MNRAS, 325, 
1288 

Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 119 

Smith, R. E., Peacock, J. A., Jenkins, A., White, S. D. M., Frenk, C. S., 
Pearce, F. R., Thomas, P. A., Efstathiou, G., & Couchman, H. M. P. 2003, 
MNRAS, 341, 1311 

Song, Y.-S. & Knox, L. 2004, Phys. Rev. D. 70, 063510 

Spergel, D. N., Be an, R., Dore, O., & et al. 2006, ApJ submitted 
jastro-ph/0603449} 



14 



RUDD, ZENTNER, & KRAVTSOV 



Suginohara, T. & Ostriker, J. P. 1998, ApJ, 507, 16 
Takada, M. & Jain, B. 2004, MNRAS, 348, 897 
Takada, M. & White, M. 2004, ApJ, 601, LI 

Tegmark, M., Eisenstein, D. J., Strauss, M. A., & et al. 2006, Phys. Rev. D, 
74, 123507 

Tegmark, M., Strauss, M. A., Blanton, M. R., & et al. 2004, Phys. Rev. D, 69, 
103501 

Tinker, J. L., Weinberg, D. H., & Warren, M. S. 2006, ApJ, 647, 737 
Tyson, A. 2005, in ASP Conf. Ser. 339: Observing Dark Energy, ed. S. C. 

Wolff & T. R. Lauer, 95 
Vale, C. & White, M. 2003, Apj, 592, 699 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, Markevitch, M., Murray, 

S. S., & Van Speybroeck, L. 2006, ApJ, 640, 691 
Wechsler, R. H., Zentner, A. R., Bullock, J. S., Kravtsov, A. V., & Allgood, 

B. 2006, ApJ, 652, 71 



White, M. 2004, Astroparticle Physics, 22, 211 
White, M. & Hu, W. 2000, ApJ, 537, 1 

Wood-Vasey, W M., M iknaitis, G., Stubbs, C. W., & et al. 2007, 

ApJ submitted ( astro-ph/0701041 1 
Zehavi, I., Weinberg, D. H., Zheng, Z., & et al. 2004, ApJ, 608, 16 
Zeldovich, Y. B., Klypin, A. A., Khlopov, M. Y., & Chechetkin, V. M. 1980, 

Soviet J. Nucl. Phys., 31, 664 
Zentner, A. R. 2006, Int. J. Mod. Phys. D, In Press, lastro-ph/0611454 1 
Zhan, H. 2006, Journal of Cosmology and Astro-Particle Physics, 8, 8 
Zhan, H. & Knox, L. 2004, ApJ, 616, L75 
Zheng, Z. 2004, ApJ, 610, 61 



